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A dark sector resembling the standard model, where the abundance of matter is explained by 
baryon and lepton asymmetries, and stable constituents bind to form atoms, is a theoretically 
appealing possibility. We show that a minimal model with a hidden SU(2) gauge symmetry broken to 
U(l), with a Dirac fermion doublet, suffices to realize this scenario. Supplemented with a dark Higgs 
doublet that gets no VEV, we readily achieve the dark matter asymmetry through leptogenesis. The 
model can simultaneously have three portals to the standard model, through the Higgs, nonabelian 
kinetic mixing, and the heavy neutrino, with interesting phenomenology for direct and collider 
searches, as well as cosmologically relevant DM self-interactions. Exotic bound states consisting of 
two fermions and a doubly-charged vector boson can exist in one phase of the theory. 


Dark matter (DM) from a hidden sector has been a 
popular alternative to supersymmetric weakly interact¬ 
ing massive particles in recent years nig. A widely 
studied example is dark atoms, where the DM consists 
of two species with opposite charges under an unbroken 
U(l)ft, hidden sector gauge symmetry [SHS]. This class 
of models presents rich possibilities for direct detection 
[iHHni, as well as cosmological imprints [TIHT^ . If the 
hidden photon has kinetic mixing with the normal pho¬ 
ton, the dark constituents acquire electric millicharges 
m, leading to further constraints and prospects for de¬ 
tection pTk25] . 

Simplified models of atomic dark matter are easy to 
construct, consisting of just two fermions and the gauge 
boson in the hidden sector, but such examples are nec¬ 
essarily incomplete descriptions of the new physics re¬ 
quired. First, it is desirable for the DM to be asymmetric, 
otherwise the long-range U(l)/j interaction would leave 
too small a relic abundance unless the DM mass exceeds 
~ 400 GeV [53]^. Simplified models do not explain the 
origin of the asymmetry. Second, the U(l);j gauge in¬ 
teraction leads to a Landau pole at high energies, so it 
would be desirable to find a more UV-complete version of 
the theory. Third, dark constituent millicharges greater 
than ^ 10 “^e (of interest for collider searches) require the 
atomic constituents to be nearly equal in mass, which is 
a rather ad hoc requirement in the simplified models. In 
this work we present a model that is still relatively simple, 
but addresses both of these issues, and makes a number 
of interesting experimental predictions. It relies upon 
breaking a nonabelian (hence asymptotically free) gauge 
symmetry SU(2)?, down to U(l)/j to explain the origin 
of the massless dark photon. The approximate equality 
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^ For lower masses the DM self-interactions violate bounds from 
structure formation. This argument assumes that the DM re¬ 
mains ionized, which turns out to be valid for the gauge cou¬ 
pling strength needed to get the right relic density from thermal 
freezeout. 


of the dark consituents, if desired, can be explained as a 
remnant of the gauge symmetry. 

There have been many proposals for mechanisms that 
link the asymmetries of the hidden and visible sectors. 
In general, they tend to be complicated. A notable ex¬ 
ception is to use the out-of-equilibrium decays of heavy 
neutrinos to generate both asymmetries via leptogenesis 
and its analog in the hidden sector [3 . We adopt 

this approach here. 

The model presents opportunities for direct detection, 
either through Higgs portal interactions or nonabelian 
gauge kinetic mixing. The latter can arise through a 
dimension-5 operator involving the triplet Higgs field 
that breaks the SU(2);i gauge symmetry |3I]. This re¬ 
sults in electric millicharges for the dark matter con¬ 
stituents, that normally must be very small to avoid di¬ 
rect detection, but can be sizable if the dark constituents 
have equal mass, which is a symmetry limit of the the¬ 
ory presented here. Moreover the self-interactions of the 
dark atoms can be of the right magnitude for addressing 
problems of small-scale structure formation in standard 
noninteracting ACDM cosmology. 

In the following we introduce the model (section 
and then estimate the dark matter and baryon asymme¬ 
tries that can arise in a generic scenario for leptogene¬ 
sis (section]^. Limits from direct searches are worked 
out in sectionIn section]^ we consider the region of 
parameter space in which the vector bosons are stable, 
leading to a markedly different dark sector. In sect.j^we 
discuss constraints pertaining to the ionization fraction, 
dark atom self-interactions, and searches for millicharged 
particles. Conclusions are given in sect. 


1. THE MODEL 

The new-physics content of the model (summarized in 
table 1^ is a hidden SU(2)?, gauge boson with field 
strength a real scalar triplet (j) that spontaneously 
breaks SU(2);j by getting a VEV, a scalar doublet 77 that 
does not get a VEV, two Weyl fermion doublets ijjf (with 
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gauge index a and flavor index i) and the heavy right- 
handed neutrinos Nj that also interact in the usual way 
with the standard model neutrinos. An even number 
of fermion doublets is required to avoid Witten’s global 
SU( 2 ) anomaly [ 32 ] ■ They can be combined into a Dirac 
doublet fermion = {iPil,4’2r) where the conjugate is 
defined as = cf2'T2{'4’2l)*^ *-e., the epsilon tensor is 
applied both to the spin and to the SU(2);i gauge indices. 
Without loss of generality the VEV of (j) can be rotated 
to the 3 rd component, ((/>“) = ( 0 , 0 , cr). 


particle 

VEV 

50,++.— 

(0, Q,a + 8) 

g+’ 


4/+’- 

0 

Spin 

1 

0 

0 

1 

2 

1 

2 

1 

2 

SU( 2 );, 

3 

3 

2 

2 

2 

1 

U(l). 

0 , +2, -2 

0 

+ 1, —1 

+ 1, —1 

+ 1, —1 

0 


TABLE I: New particle content in the model, showing the 
Lorentz, hidden SU( 2 ) and hidden U(l) (after breaking of 
SU( 2 )h —>'U(l);i) quantum numbers. 


The relevant terms in the Lagrangian are 

C = - + (1) 

-I- ^{ 0 {yi+iy 2 lh){$-T)'l> ( 2 ) 

- \D^ 0 -V{H 00 

- 0L'n)y^PRNj+Yc. 

where the covariant derivative is D 0 °' = 80 °^ — 

g Cabc or D 0 = [ 8 ^ - i{gl 2 )B^ ■ rjT', g is the 

SU( 2 )^ gauge coupling, and is the Standard Model 
hypercharge gauge field strength. In the Yukawa interac¬ 
tions with the sterile neutrino we use the Weyl fermion 
notation since the analogy to leptogenesis via neutrino 
physics is more clear in this way. 

The triplet scalar VEV breaks SU( 2 );i to U(l)?i, me¬ 
diated by the massless gauge boson iJg, while B^ = 
[B^ ± iB'^)l\f2 obtain mass tob = g<J. The upper 
and lower components ^'1^2 of the fermion doublet are 
also charged under the U(l)/j (with half the charge of 
B^). Their masses are split by the Yukawa interaction, 
mi,2 = {{m^ ± 2/1 cr)^ -I- (2/20')^)^/^. We used the freedom 
to perform a chiral rotation on 'k so that is real (has 
no 75 component). 

In section | 2 . 4 | we discuss the decay of the scalars 
through r] —>■ ijji' via the dimension -5 operator 

0 LV {H^LD + h.c. ( 3 ) 

where Mj is the heavy neutrino mass (in a basis where 
its mass matrix is diagonal), y^ is the neutrino Yukawa 
matrix, H is the SM Higgs doublet, and Lk are the lepton 
doublets. 

We will initially consider the case where decays T'l —>■ 
H++ 4'2 are not kinematically allowed. They would lead 
to a dark sector consisting of stable vector bosons 


0 “ 


I 



FIG. 1 : Loop contribution to the nonabelian kinetic mixing 
operator. 


and fermions. (The alternative case in which these 
decays are allowed is considered in section]^) This leaves 
two species of stable dark matter, the Dirac fermions 
^1 = {'<P\l02r)'^ and ^-2 = 01l02rV with charges 
±1 under the unbroken U(I)/j. The long-range force me¬ 
diated by the dark photon B3 = 7' causes the symmetric 
component of the DM densities to be at least partially 
depleted by annihilations, and the asymmetric compo¬ 
nents of d'l 2 to bind into dark atoms. The efficiency of 
these processes depends upon the gauge coupling g and 
the dark atom mass toh, as we will discuss in section]^ 

For simplicity we impose a softly broken U(I) symme¬ 
try under which ipi —)■ g —)■ which forbids the 

interactions (jp^Nj, with fj = T2g*■ The symmetry is 
broken by the Dirac mass term, which takes the form 

- +^2flV’iL)+h.c. ( 4 ) 

If the symmetry were exact, then the subsequent decays 
g ^ ip mediated by Ni would completely erase any pro¬ 
duced DM asymmetry. However the chirality flips in¬ 
duced by the mass term prevent this erasure, as we will 
explain in more detail in section]^ There is an unbroken 
discrete Z2 remnant of this symmetry, where ipi —^ —ipi 
and g —>■ —g, that ensures the stability of the dark mat¬ 
ter. 

The potential V is assumed to give rise to the VEV of (p 
and it generically also includes the Higgs portal coupling 
^Xh^\H\'^(p^■ Once (p gets its VEV, the nonabelian kinetic 
mixing operator can be written as 

-isine'B^^Y'^^' (5) 

where sine = 2 cr/A. It could arise from integrating out a 
heavy vector-like fermion x that carries hypercharge and 
transforms as a doublet under SU( 2 );,. The interaction 
yxX^Pa^^aX leads to the diagram in flg.[^ implying A“^ ^ 
ggiyx/mx where gi is the hypercharge coupling. The 
kinetic mixing gives rise to electric millicharges ±617 = 
±ee for the fermions 'I'l 2- This or alternatively the Higgs 
portal interaction allows for direct detection of the dark 
atoms, as we discuss in section]^ 
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2. ORIGIN OF DARK MATTER ASYMMETRY 


Our setup allows for heavy neutrinos to decay in a CP- 
violating manner into an excess of dark matter versus its 
antiparticles in close analogy to leptogenesis. The struc¬ 
ture of Yukawa couplings is similar to that of neutrinos 
except that we have only two light fermionic DM species, 
as compared to the three light neutrinos. The dark 
Higgs doublet 77 does not have a VEV, so it also gets an 
asymmetry, which will be determined by those in 'I'i. 

The asymmetry in the decay of the jth heavy neutrino 
into 7/1*77 versus ^^77* (recall that t/i^ denotes the Weyl 
doublet states) is given by 


el;, = 


r(Y, ^ 7 /;* 77 ) - T{N, ^ 7 /-, 77 *) 


T{Ni -)■ any) 


( 6 ) 


= -T 

Stt 




Im 




{ylvi^ + ylyu)n 


Im 


{ylyu)k]y'‘^y^* 


{y\py^ + ylyv)n 






where g{x) = y/x[l/ {1 — x) + 1 — {1 + + 1 /x)] 

and g'{x) = ^/x/{l — x). This differs from the standard 
leptogenesis expression because the denominator must 
take into account decays of Nj both into neutrinos and 
dark matter, and there is a mixed term of order 7/^ 
from the self-energy correction of Nj by the SM Yukawa 
interaction. 

For definiteness, we will focus on decay of the lightest 
heavy neutrino Ni. In the simplest scenario of leptoge¬ 
nesis, where Mi <C M2,3 and the reheat temperature is 
in between, Mi < Trh < M2,3, this is the only relevant 
decay since the heavier neutrinos are not present. In 
this case the functions in eq. ® can be approximated as 
g = — 3 / 2 y^ and g' = —I/y^icwith x = (M2/Mi)^ ;g> 1. 

Initially, we can expect independent asymmetries Yi,2 
for 7 />i and 7/12, where = {n^. — n^^js is the dark 
matter to entropy ratio, since However the 

Dirac mass term takes the form ipi a2T2'tjj2, which implies 
that mass effects will cause the asymmetries of ifji 
'ip2 to become equal and opposite. This projects the net 
asymmetry of the fermions onto the difference between 
the initial ones, = Yi — Y2, at temperatures where 
the helicity-flipping interactions due to come into 
equilibrium. 

On the other hand, the 77 boson gets a different asym¬ 
metry, proportional to + e]p . Eventually it will decay 
into ij^i. For simplicity, we consider the case ~ 

Then not only does the initial asymmetry in 77 tend to be 
small, but so also is its contribution to the final asym¬ 
metry in ipi, and we can estimate the net asymmetry in 
7/) from iVi decays as 

e^i ^ (7) 

The sign difference is in contrast to the CP asymmetries 


for decays into neutrinos, e^i = e/,* familiar from lep¬ 

togenesis. 


2.1. Dark matter asymmetry estimate 


The initial asymmetries depend upon an efficiency fac¬ 
tor Kj, that quantifies the amount of washout (see for 
example | 33 j for a review). The contribution from Ni 
decay is 


Yj, = 


45 (L'lji Kj} 

TT^ 


( 8 ) 


where k,^ = min( 0.25 (tti^/tti^i)^'^, 1 ) with 771^1 = 
2(77j,7/^)ii7;^/Mi, 777* = 10 “^eV and v = 174 GeV. 
The Higgs VEV v has no direct physical relevance for 
the dark matter abundance, but 777^1/777* gives r(Vi —>■ 
ipr]^*'>)/H{Mi) (the ratio of the partial decay width to 
the Hubble rate), just like rhu^/m^ does for the decays 
into vh. The dark sphalerons associated to the SU( 2 )^ 
gauge interactions have the same effect as (an increase in) 
the Dirac mass term for tk and therefore do not require 
additional consideration for the dark asymmetry. 

If there is no hierarchical structure to the couplings 7/^ 

and their phases are large, we can estimate {ylpytp)kj ~ 
{ylpy,p)kj or its imaginary part by some average value 
Further defining yl = {ylyij)ii and assuming that 
the terms of order y^ in the numerator of ([^ can be 
estimated as yl, we find that the CP asymmetry for 7/; is 
of order 


yif’ f ^ 4~ 2^ A 

Stt ^/x \ 1 -I- r y 


( 9 ) 


where we define r = y^/yl, and assume that ^ —e]p 
in (|^. It is evident that eq. Q has only mild depen¬ 
dence upon r. Combining with the efficiency factor Kj; 
(where we approximate the exponent 1.1 by 1) leads to 
the estimate^ 


V/, ^ 1.4 X 10-^2 


Ml \ 

IQio GeV J 



( 10 ) 


ignoring r dependence. 


2.2. Baryon asymmetry estimate 

We wish to explain the baryon asymmetry simultane¬ 
ously with that of dark matter. Analogously to ([^, it is 


^ This is val id f or parameters such that Kj, < 1 hence m > 4m*. 
Using eq. |l3| this implies > 10~^(Mi/10^'^ GeV). We will 
assume that this restriction holds in the following. 

















4 


given by 


79 


45 

7r4 


( 11 ) 


where the prefactor 28/79 is due to redistribution of the 
initial lepton asymmetry into baryons via sphaleron in¬ 
teractions. The CP asymmetry is defined as = 
in the usual way for leptogenesis. Similarly to 
our estimate in we expect the well-known D-I bound 
[Hi] to be modified by a function of r = 




3 Ml Ami 


IGtt 

^ 10 "® 


TO ^3 

Ml 


IQio GeV 


i + i^ 

1 -1- r 

= 10 "® M 


10 


( 12 ) 


where Aml^ = — m^^, which we assume to be 

= = (0.05 eV)^. Again the dependence upon r is 

mild, and we will ignore the effect of the DM Yukawa 
coupling on leptogenesis in the visible sector. To esti¬ 
mate the efficiency factor K,y = 0.25 (m^/rh^), with = 
iyly<^) iiu^/Mi, we use the Casas-Ibarra parametrization 
of yu, 


iylyi^)ii 



Rkjml(^ C/Ji 


10 "® Mio 


(13) 


(the same result as eq. (lil) where U is the PMNS 
matrix and R is an arbitrary SU(3) transformation. 
We assumed that R\f.MkRkj ~ Mi since we take the 
heavy neutrino masses to be of the same order, and 
|tli 2 pW!V 2 + |lli3pnr!^3 = 0.003 eV (taking to be much 

less than the solar neutrino mass splitting). This gives 
Kj, = 1/12 and 


Yb ^ 1.4 X 10-1° A/io (14) 


where we have introduced a parameter e^i to quantify 
how much falls below the D-I bound, be., Cm is |ey,i| 
over its maximum value. Equating Yb to its measured 
value, we find em Mio = 0.7. 


is interesting for direct detection and consistent with our 
prejudice for the new physics scale to not be far below the 
weak scale. It is interesting that the same mass scale is 
also consistent with the observed baryon asymmetry for 
generic choices of the neutrino CP asymmetry, eoi < 1. 

Notably absent from our estimates is any explicit de¬ 
pendence upon the Yukawa couplings and This 
is because of the cancellation between the CP asymme¬ 
try e and the efficiency factor k, which only occurs for 
couplings such that k < I. We verified this condition for 
Ki, = 0.08. It is also satisfied by so long as > yijYl. 
We will make this technical assumption to keep the anal¬ 
ysis simple. For smaller values of there would be a 
suppression of Y^ and the need for correspondingly larger 
values of toh- 

As an example, we take = 0.65, y'^ = y1 = 10"®, 
ruH = 83 GeV, Mi = lO^® GeV, M 2 = y^Mi = 2 x 
10^^ GeV. Larger or smaller values of mn can be obtained 
by adjusting M 2 /M/, using eq. (15). 


2.4. Decay of dark scalars 


An interesting feature of our model is that the seesaw 
mechanism produces the new dimension-5 operator ([^ 
that allows the dark scalars y to decay [29]. When the 
SM Higgs takes its vacuum expectation value, this allows 
the rj to decay directly into vil}. The decay rate is of order 


r 


ylylm^v'^ 

87rM2 

3 X 10"° s"^ 


■( 


150 GeV/ 


(16) 

(17) 


where for the numerical estimate of P we used the exem¬ 
plary values specified at the end of the previous section 
(ignoring the mass of tp in the phase space integral). 

Such decays must occur sufficiently early so that the 
decay products are fully thermalized before they can dis¬ 
tort the GMB. Ref. [3S] shows that this occurs if the life¬ 
time is below ~ lO^^s. Eq. (16) implies that our model 
easily satisfies this bound. 


2.3. DM to baryon constraint 


3. DIRECT DETECTION 


We can combine the above results to get a con¬ 
straint on the model parameters from the known ratio 
of baryon and dark matter energy densities, Q.b/^dm = 
knpYB/{miYY^) = 0.18. Here tor = + ^2, the mass 

of the dark atom (neglecting its binding energy). Then 
we find m-glmp = 166 Cm (a;/10)^/^. We can eliminate 
Edi using eq. (14) and Mi using © to obtain 


Toh 


= 560 Er 


p.1/2 

10 


360 

w~ 


10 


j.1/2 

IcT 


(15) 


Eq. (15) reveals part of our motivation for the choice 
Ml ^ 10^® GeV: it gives dark atom masses in a range that 


There are two portals through which our dark atoms 
to interact with nuclei. The kinetic mixing allows for 
photon exchange, which has been discussed in refs. |211 
[HH] . The ensuing constraints on the electric millicharge e 
are weakened for atoms compared to ions because of the 
screening of electric charge. In the special case where 
mi = m 2 this screening is perfect, and the interaction 
becomes magnetic dipole, further weakening the limits 
m- Here we extend results of ref. m for the mi = m 2 
case to higher DM masses. 

In addition, there is the Higgs portal induced by mix¬ 
ing of h and (ps through the operator ^Xh(i,\H\‘^\(p\'^. (p^ 
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FIG. 2: Left: CRESST-II, CDMSlite and LUX limits on millicharge e of dark atom constituents, with constituent mass ratios 
m 2 /mi = i? = 2, 4, ■ • ■ , 10 as indicated, for photon-mediated scattering of dark atoms on protons. For clarity, only the most 
constraining limit is shown for any DM atom mass mu- Gauge coupling is set to Og = aion, eq. (191 for solid curves, and fixed 
at Qg = 0.06 for light dashed curves. Right: Corresponding limits on F = yi 6\1 — m^/m^|F’,/,(0) from Higgs portal scattering, 
where 9 is the c^-Higgs mixing angle. 




FIG. 3: Direct detection constraints on kinetic mixing parameter e versus dark atom mass rriH for case of equal-mass constituents 
mi = m 2 = mu/‘i, when interaction is magnetic inelastic. Hyperfine mass splitting is chosen as a function of mu as described 
in text. Solid and dashed curves refer to choice of a.g as in fig. Left: low mass region; right: larger mass region. 


interacts with the dark atom constituents through the 
operator ^(yi+fy 275 )(<?’^)^ that splits the 'I'i _2 masses. 


3.1. Kinetic mixing portal 

3.1.1. Unequal-mass constituents 


As discussed in ref. m, the mutual screening of the 
electric charges of the dfi and 'I '2 constituents results in 
a scattering matrix element where the 1 /q^ of the photon 
propagator is canceled by q'^ in the form factor for the 
charge density. The cross section for scattering of dark 


atoms on a proton is 


= Itt 


a;2 £2 ^2 


q ;2 £2 u 2 

= An^^foiR) 


4 4 

ag m 


(18) 


H 


where fj-n is the reduced mass of the dark atom and nu¬ 
cleon system, and fo{R) = (1 + — 1)^. Here 

we have generalized the result of ref. [21j where the ap¬ 
proximation of large R was made. The expression (18) 


is valid if R is not too close to 1. The question of “how 
close?” is discussed below. 

For i? ^ 1, the resulting upper limits on e is illustrated 
in fig. ^a) showing the most constraining limit from the 
LUX 1^, CRESST-II [57] or CDMSlite [55] experiments, 
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at any given dark atom mass mn • In section we will 
see that the requirement of sufficiently small ionization 
fraction in the dark sector leads to the constraint 


More quantitatively, the actual cross section for a given 
scattering event is Cp = ap^F^q, v) and the detection rate 
is proportional to 


as > aion = 5 X (19) 

where 


f2(E) =E + 2 + E-^ ( 20 ) 

The solid curves are derived for the parameter choice 
which saturates this bound, ag = aion, and F ranging 
from 2 to 10, while the dashed ones assume a fixed value 
of ttg = 0.06. This value satisfies the constraint Og > aion 
over the entire range of E and mu shown on the plots. 
(The unusual sensitivity of the solid curves to light DM 
masses is due to the decrease of ag = aion with toh, and 
consequent increase in the dark Bohr radius, leading to 
larger cross sections.) The nominal constraints from the 
experiments are weakened by factors of (A/Z)^ = 5.9, 5.2 
and 4 respectively to account for the coupling to protons 
only. For CRESST this corresponds to collisions with 
the oxygen atoms that dominate the senstivity to low- 
mass dark matter. The strongest constraints occur for 
mn = 1 — lOGeV, in the range e < 10“^° — 10“®. For 
conventional abelian kinetic mixing, such small values of 
e could be difficult to achieve since the loop diagram that 
generates it is not surpressed by any large mass scales, 
since in this case the kinetic mixing operator is marginal. 
However for nonabelian kinetic mixing, e is suppressed by 
the mass of the heavy particle in the loop, as well as 
its Yukawa coupling ^ . For example if the couplings 
described below eq. ^ are = 0.1, gi = g, E = 10, 
ctg = aion and cr = 30 GeV, we require rriy, > 3 x 10^^ 
GeV to satisfy the LUX bound on 10 GeV dark atoms. 


3.1.2. Equal-mass constituents 

For i? = 1, there is perfect screening of charge because 
of the complete overlap of the wave functions of the two 
constituents, and the magnetic dipole interaction that we 
have neglected in (181 becomes important. This case was 
considered in detail in ref. [5T] . The magnetic scattering 
is inelastic because of the hyperfine transition of the dark 
atom, requiring energy 5E = mu, hence a minimum 

DM velocity of Umin = 9/(2^w) SE/q > ^2 SE/qiN 
for momentum transfer q and dark atom-nucleus reduced 
mass gN- There is a q- and u-dependent form factor F = 
( 90 / 9 )^ ~ of order unity for typical 

values u ~ uo and q qo, as along as vq > Umin- The 
cross section on protons is of order 


^p,0 = 




( 21 ) 


in that case, where is the proton-atom reduced mass. 


/■E„ 


E^Z'^ I cIEr - f{v)ap 

.In) . „ V 


(X Z^apfilp, 


( 22 ) 

(23) 


with 

P _ 1 2 

Eniin — 

Einax = Pmax/(2mjv) 

Pmax = \/ gliVesc + Vq)"^ - 2 6EgN -|- flN {Vesc + fo) 

f{v) cx (24) 

/-WeBC J3- 

If= cIEr ^/(F)F(g,u). (25) 


Here Ve is the Earth’s speed relative to the DM halo, vq « 
220kms“^ is the mean DM velocity, Uesc ~ 450kms“^ 
is the approximate escape velocity of the DM halo (we 
see no significant variation in the results for values in the 
range 400 — 500 kms“^), and Egsc is the maximum recoil 
energy from a DM particle with the escape velocity. 

We compare the rate for our model to that of generic 
DM scattering with a constant cross section cr„, for which 
the corresponding expressions are 


E oc A^ct„/o, 




In = 


dEf 


d^i 


fiv)- 


(26) 

(27) 


where Emix = Emax evaluated at 6E = 0. Therefore 
the magnetic inelastic cross section (21) is bounded from 
above as 


^P. 0 < 


lim 


(28) 


where (TuXim is the experimental upper limit on the cross 
section for a generic DM model. Notice that the arbitrary 
quantity (vo/qo)'^ appears in the same way on both sides 
of (28) and hence can be divided out. 


Although the gauge coupling ag does not appear in 
(1^, the mass splitting SE = agmu/Q depends upon it. 
For definiteness, we have chosen the value ag = aion in 


eq. (19) from the requirement of sufficiently small dark 


ionization fraction. This fixes SE as a function of mu. 

We plot the ensuing limits on e in figure]^ using results 
from the LUX [35], SuperCDMS [33], CRESST-H [37] . 
and CDMSlite [33] experiments. The mass dependence 
of ag = Q!ion changes the shape of the exclusion curves 
relative to those on the cross section itself, since the mass 
splitting SE rises rapidly with mu, nullifying the signal 
for mu ^ 100 GeV. 
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FIG. 4: Ratio of the magnetic inelastic and elastic cross sec¬ 
tion for scattering of dark atoms on protons as a function of 
the constituent mass ratio R (its deviation from unity), for 
dark atom mass thh = 10 GeV and mass splitting described 
in section [3.1.2l 


3.1.3. Transition from R = 1 to R> 1 

We have noted that inelastic magnetic transitions dom¬ 
inate for equal constituent masses, R = 1, while elastic 
charge-charge interactions dominate when i? > 1. One 
may wonder how sharp the transition is between the two 
regimes; how small must i? — 1 be for inelastic transi¬ 
tions to dominate? We have calculated the ratio of the 
two cross sections as a function of R for a particular 
value of toh = 10 GeV as an example, taking the mass 
splitting 6E as described above. The result is graphed 
in fig. 1^ which shows that only for R < 1.0001 do the 
inelastic transitions dominate. Hence the most natural 
situation corresponding to this case is that where R = 1 
exactly. There are two limits in our model that give 
R = 1: either yi = 0, or = 0. The latter is a point 
of enhanced SU(2) flavor symmetry for the two chiral 
(doublet) fermions. 


3.2. Higgs portal 


The interaction of dark atoms with the Higgs through 
(p 3 -H mixing also undergoes screening because of the cou¬ 
pling Ta which has opposite sign for 'I'l and '^ 2 - At low 
energies, the dark atoms can be described by a Dirac 
field H whose coupling to the virtual (^3 or h carrying 
momentum q is given by the amplitude 

yiu^u^ F^{q) (29) 

We have neglected the 2/2 contribution that is suppressed 
by the dark matter velocity. By matching onto the scat¬ 
tering amplitudes in the high-energy theory, we infer that 


= - 

mn 


m 2 


mi 


_{l + \q^alY 


( 30 ) 


with Gi = (agmi)~^. Thus the coupling vanishes in the 
limit i? = 1 (toi = m 2 ). If 0 is the mixing an¬ 
gle, then the amplitude for scattering of dark atoms on 
nucleons is 


M = 2/1 Uh(P3)Mh(Pi) ■ Un(P4)Un(P2)(^^^^^) (31) 


X 




F^{q) 


where (ynmnlv) with 2/n — 0.3 [10] is the coupling of the 
Higgs to nucleons. 

If ag is not too small, we can take the q — 0 limit of 
the form factor. In this case the cross section for dark 
atom-nucleon scattering is 

= ^^2 [yi ynmjiqinM ^ F.^p(0)] {m^ 771 ^ ) ( 32 ) 

in terms of the H-nucleon reduced mass, and taking 
0^1. The LUX upper limit on the dimensionless com¬ 
bination F = yi 9\1 — m^/m^|F. 0 (O) is plotted in fig. 
The Yukawa coupling 2/1 is related to the mass splitting 
in the dark sector since m\ — m\ = 2/1 Moreover it 
is straightforward to show that F^{Q) = {m\ — 

If 7710 < mh then F = . We expect 

7770 ^ (T, similarly to mh ^ u in the visible sector, and 
6 < 0.01 to satisfy LEP constraints [H] on mixing of a 
light scalar with the Higgs. The largest dark atom mass 
range for saturating the LUX bound shown in fig. with 
I 2 / 1 I < 1 is 777 h ^ 70 GeV. 


4. STABLE VECTOR BOSONS 


Up to now we have assumed that the 'I'i-'I '2 mass split¬ 
ting is sufficiently small to prohibit the decay 41^ —>• 
B 4')'', corresponding to the condition 


\yi\ < 


7771 + m 2 
47770 
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(33) 


However this need not be the case, and the model is also 
compatible with a universe where charge neutrality in 
the dark sector is acheived by having two 41)^ particles 
for every B . This leads to a very different kind of dark 
atom that is reminiscent of the H 2 molecule, except that 
the two “protons” are bound together by a single charge 
—2 “electron”. We will refer to these variant dark atoms 
as H 2 . In the absence of fine-tuning, the stable vector 
boson is typically lighter than 471 , prompting us to define 
the ratio 


i?2 = — > 1 

777b 


(34) 


in analogy to R = m 2 lmi for H atoms.^ 


^ To get the opposite situation where rriQ > mi, we need (p'cr)^ > 
{rrifjj — yicr)^ + (y 2 cr)^- This requires not only y 2 to be small, but 
also an accidental cancellation between m.^ and yicr. 















FIG. 5: Left: energy obtained from variational method as a fnnction of R 2 for the trial wavefunctions for H 2 bound states 
tpH 2,3 and tpH 2 ,p 2 - Right: Corresponding values of the parameters a, b (shown in dimensionless combinations with agms) that 
determine the spatial distributions of the wave functions, for the s-wave. 


4.1. Bound States 

To verify the existence of the 3-body H 2 bound states, 
we make some ansatze for its wave function and use the 
variational method to prove that the energy is minimized 
at a negative value. We consider trial wave functions 
where the positions of the three particles are given by 

= ±A/2, xb = r (35) 


'kH 2 ,pi is always more weakly bound than the other two. 
Moreover the s-wave has lower energy than p2 only for 
^2 ^ 40; for i ?2 > 40 the p2 state is lower, as shown in 
fig-d Taking as an example the values mi = 60 GeV, 
i?2 = 10, Og = 3 X 10“^, the three-constituent atoms 
have binding energies of approximately E « — 15MeV. 


4.2. Direct Detection 


i.e., we work in the center-of-mass frame of the two iki 
particles, with A being their relative separation. In anal¬ 
ogy to the H 2 molecule, it could be expected that the 
wavefunction for A is approximately that of a 3D har- 
monic oscillator, e~ ' for some scale b. For simplicity 
we take the wave function for r to be hydrogen-like, 
for some other scale a. We consider three possible states, 
an s-wave and two p-waves, 

V^H2 .s(gA) = (36) 

V^H2.pi(GA)=iVpir, 

V^H2.p2(r, A) = Ap2 A, 


where (Aj,) is the ^-component of r (A). 

It is convenient to work in the analog of atomic units 
by rescaling to dimensionless coordinates r = r'/(ogTOs), 
A = is! j(oLgmB)- Then the Hamiltonian can be written 
as H = {agmB)H', where the dimensionless H' is 




^ IF' 


r'± AV2| 


(37) 


Dark H 2 atoms interact similarly with nucleons rela¬ 
tive to our treatment for H atoms in section but there 
are some qualitative differences, due to the more compli¬ 
cated wave function. In particular, there is no longer any 
special case like i? = 1 for H atoms in which the electric 
millicharge clouds of the constituents give exactly cancel¬ 
ing contributions to the total charge density. This can be 
seen by computing the form factor, which is the Fourier 
transform of the charge density 


p{x) = / d'^A |d>(r, A)| 


X S(x ± A/2) — 2S(x — r)^ 


Using V'H 2 ,s from eq. (36), the form factor is 

^(®) = 2 + (l + aV/4)2) 

= ^ (l6"“ 


(38) 


(39) 


By minimizing the expectation values E = 
(■ 0 H 2 I^I' 0 H 2 ) with respect to a,b and varying over 
a range of i ?2 values, we find that bound states (having 
E < 0) exist for all three trial wave functions, but 


where the approximation is for low momentum transfer 

q- 

In computing the cross section for scattering on pro¬ 
tons, the factor of in the form factor cancels the l/q"^ 
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FIG. 6 : Direct detection constraints on kinetic mixing pa¬ 
rameter as in fig. but for H 2 atoms with mi/ms = R 2 = 


2,4,. •• ,10. 


of the propagator like before, giving 


CTp = 




(40) 


at low momentum transfer. (The normalization can be 
deduced by considering the limits a = 0 , & —>■ 00 or vice 
versa where the usual Feynman rules for the amplitude 
with no form factor apply.) The direct detection lim¬ 
its from LUX [55], CRESST-II [57] and CDMSlite [55] 


through the kinetic mixing portal are shown for various 
values of R 2 in figure assuming ag saturates the con¬ 
straint (51) from ionization of H 2 atoms that we will 
derive in the next section. (We also show the constraints 
for the fixed value of Og = 0.01 as dashed curves.) Unlike 
with the dark atoms, the form factor never vanishes for 
any value of R 2 (since b is always < 4a). 

For the Higgs portal, we follow the procedure in sec¬ 
tion!^ The amplitude and cross section are 


M = WH2(P3)mH2(15i) • Un{p4)u„{p2) ^ ^ 

^ (41) 

\ '"'h / 

1 2 
<Jn = —y [{yiF^{0) + g Fb{ 0)) ynrrinPnd] 

TTV 

(42) 


We have again made the approximation 0 <C 1 and as¬ 
sumed a small momentum transfer, pn is the H 2 -nucleon 
reduced mass, and ?/„ — 0.3 is the Higgs coupling to nu¬ 
cleons (modulo rrinlv). The form factors are given by 


F^{q) 

Fsiq) 


2r^^_{,2^2/32 

TOH 2 

ruB 1 


2toh 2 (1 + ia2g2)2 


(43) 


Redefining F = {yiF^{0) + g FBi0))0\l - the 

constraint on F from the LUX, CRESST-H and CDM¬ 
Slite experiments takes the same form as was previously 
shown shown in figure]^ (right), where tor is reinter¬ 
preted as TOH 2 ■ 


4.3. Neutron Star Constraints 


Tight constraints exist on the cross section for asym¬ 
metric bosonic dark matter scattering on nucleons from 
the existence of long-lived neutron stars [42l [45] . If the 
rate of dark matter accretion is large enough, it can col¬ 
lapse to form a black hole that would consume the pro¬ 
genitor, on time scales shorter than the ages of neutron 
stars observed in globular clusters. In our model it is 
important that we have only one kind of stable bosonic 
dark matter consituent carrying dark U(l);i charge. In 
the case of H atoms with only fermionic constituents, the 
would-be scalar constituents decayed early in the cos¬ 
mological history, leaving no asymmetric scalars. For 
H 2 atoms, on the other hand, the vector bosons are 
mostly bound inside of atoms that resist collapse be¬ 
cause of the degeneracy pressure of their fermionic con¬ 
stituents. The ionized fraction also resists collapse be¬ 
cause of dark Coulomb repulsion. In contrast, in a model 
containing two species of bosons carrying different U(l)?i 
charges, nothing would prevent the collapse of the com¬ 
bined bosonic fluid. 

In more detail, we first note that the dark atoms re¬ 
main bound once they start to accumulate in the neutron 
star. From figure [^ the binding energy is given by 


Ub « 2 Og TUB 


2 al mH 2 
l + 2i?2 


(44) 


Using the dark ionization constraint (51), we find that 
Eb > 130 eV even for the extreme parameter choices 
mH 2 = 1 GeV, i ?2 = 100, which is higher than the 
temperature of the star, of order 100 eV [33]. Moreover 
fermions within a neutron star are supported by their 
degeneracy pressure, given by 


P = 


(Stt^ 

5 m .0 


5/3 


(45) 


where n is the number density. A larger fermion mass 
decreases the pressure, and therefore the dark atoms will 
tend to remain bound. 

As for any ionized bosons that accumulate within the 
neutron star, their repulsive self-interaction greatly weak¬ 
ens the bounds on scattering with nucleons by prevent¬ 
ing their collapse into a black hole. Ref. [3S] finds that 
a repulsive scattering cross section exceeding 10 “®° cm^ 
is sufficient to avoid neutron star constraints for ms < 1 
TeV. In our case the cross section corresponding to dark 
Rutherford scattering is infrared divergent, but if we 
make it finite by multiplying da/dFl by (1 — cos 0 )^ (thus 
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taking into account only scatterings with significant mo¬ 
mentum transfer), it is of order a^g/m\ > 10“^^ cm^, 
where we used ( [^ and ms < 100 GeV. This satisfies 
the requirements of [45j by many orders of magnitude. 


5. OTHER CONSTRAINTS 

Dark atoms, dark matter with millicharges, and mod¬ 
els with asymmetric dark bosons are subject to further 
constraints from cosmological, astrophysical and labora¬ 
tory probes. Here we discuss those coming from dark 
recombination, self-interactions of the dark matter and 
accumulation in neutron stars, and searches for mil- 
licharged particles. 


5.1. Dark ions 


If the constituents of the hidden sector fail to combine 
into atoms, they can scatter very strongly with each other 
through the dark Coulomb interaction, contradicting the 
normally assumed properties of collisionless cold dark 
matter. From fitting to results of ref. [4], one finds that 
the ionization fraction can be estimated as [HIIIIIIS] 


Ae ^ ( 1 +10^72 


GeV^ 




-1 


(46) 


where f 2 {R) = R + 2-\-\/R (introduced in eq. (20)), and 
^ is the ratio of dark sector to SM sector temperatures. 

In [3] it was argned that observations of the Bullet 
Cluster rule out > O.I, leading to the conservative 


lower limit ag > aion (|I9|) that we already incorporated 


in our analysis of direct detection constraints. Ref. M 
estimates that there is a factor of 10 uncertainty in (46). 


We note that this leads to only a factor of 1.8 uncertainty 
in the expression for ajon- 

The ratio between temperatures can be found using 
the relation m 




9*s,sm 9tF,D 
9%,d ^ 


1/3 


(47) 


with gtS,SM and gtS,D denoting the number of degrees 
of freedom in the visible and dark sectors, and the su¬ 
perscripts 0, dec indicating the respective values today 
and at the time the two sectors decouple kinetically. The 
temperature at which this decoupling occurs is therefore 
relevant. We find that mixed Compton scattering with 
one dark and one SM photons is the most important pro¬ 
cess for maintaining kinetic equilibrium. It goes out of 
equilibrium when H = (av), leading to the estimate 


1.66 0 * 


T2 


g*T^ 


, Stt 


Thus mixed Compton scattering keeps the two sectors at 
the same temperature until 


Tdec 


3 X 10 ® eV / mu \ 

^2 lew/ 


(49) 


The lowest value of Tdec is obtained by saturating the 
direct detection limits on e as a function of toh, as shown 
in figs. |2][^ In the case of i? = 1 (equal mass dark atom 
constituents), this can be much lower than the dark re¬ 
combination temperature Tree, so that in fact Tdec = Tree, 
since Compton scattering is no longer efficient on neutral 
atoms. For i? > 1 on the other hand, the constraints on 
e are sufficiently strong that the decoupling temperatnre 
is limited to Tdec > 300 TeV. 

As long as Tdec ^ 1 TeV, all particles are relativistic 
except for the heavy neutrinos. We therefore use the val¬ 
ues 5 °s,sm = 3.94 [H], 5?|,SM = 106.75, g% jj = 2, and 
D = 1®- The resulting temperature ratio is ^ « 0.71. 
At the other extreme, decoupling occurs after electrons 
have frozen out. This corresponds to gfs^sM — 7-25, 


9?S^dm = 2, and ^ « 0.81. Even at the two extremes, 
therefore, the difference is minimal, and is further mit¬ 
igated by the fact that ^ is raised to the 1/4 power in 
calculating aion- We therefore adopt the value ^ = 0.71 
in eq. (191 so that aion remains a reasonable lower limit 
for ag. 

There are certain cases that can lead to a lower tem¬ 
perature ratio, with the smallest being that in which all 
dark content apart from the dark photon has frozen out 
prior to the freeze-out of the top quark, with decoupling 
occurring some time between these; in this case the dark 
temperature could be as low as 0.3. These cases, however, 
are unrepresentative and only apply to a narrow range of 
values of e. Even in the extreme case of « 0.3, the 
estimate on aion would only differ by a factor of « 0.8, 
which is smaller than the error due to the uncertainty in 
the ionization fraction. 


5.1.1. H 2 ionization 

For the case where ^'2 can decay to and the vector 
boson H, to make a rough estimate of the ionization frac¬ 
tion, we assume that recombination will typically happen 
in two steps: in the first, unbound Ti’s combine with the 
free H’s to make a 'k-H ion, while in the second these 
ions bind with a second 'ki. The first step is similar 
to hydrogen atom recombination with the substitution 
Ug —>■ 2ag due to B having charge 2. In the second step, 
the potential at long range is like that for hydrogen atom 


recombination. 

Equation (46) then becomes 

Ael = 

fl + r'l6x 10i°a4 ^ 


\ ^mirnsj 

Xe2 — 

\ '^mi[mB + mi) 


Xel + Xe2 — Xe2 


mpi 


3 m-H ’ 


(48) 


(50) 
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The constraint on the ionization fraction {Xe^tot ^ 0.1) 
from |4] is therefore 

a, > 4 X 10-3 (g|) (51) 

where / 3 (i? 2 ) = (^2 + 1/2)^/ iR 2 + R^)- 

5.2. Self interactions 


boson part of the wave function, rather than the charac¬ 
teristic distance b between the fermions, even though b ^ 
2a. This is because the expectation values are (r) = 1.5 a, 
(A/s) = 0.46. Therefore in parallel to the H atom case, 
we can estimate the elastic cross section for atom-atom 
scattering as ct = 100 = 100 a~^m^^R 2 f^{R 2 ), where 

UiR) = l + {2R2)-\ 

The gauge coupling corresponding to the desired scat¬ 
tering cross section of a/mu 2 =1.1 b/GeV is therefore 


Although standard cold dark matter is considered to 
be noninteracting with itself, there has been interest in 
variant theories where dark matter has an elastic self¬ 
scattering cross section of order lb per GeV of DM mass. 
This has been motivated by persistent discrepancies be¬ 
tween predictions of A-body simulations and observed 
properties of dark matter halos. While simulations tend 
to predict cuspy density proHles for galaxies, there is 
some observational evidence for cored profiles, especially 
in dwarf spheroidals. Simulations also tend to predict too 
many high-mass satellite galaxies accompanying Milky- 
Way like progenitors compared to observations. For a 
review of these problems and their possible resolutions, 
see ref. |4t)] . A number of studies have been done in¬ 
dicating that the small-scale structure problems can be 
alleviated by invoking dark matter elastic scattering with 
a jm ~ Ib/GeV. Dark atoms can naturally accommodate 
such large cross sections since they can have a significant 
geometric size. 

The elastic scattering of dark atoms on each other 
has been studied very quantitatively, thanks to the fact 
that the problem can be mapped onto that of normal 
atom scattering with appropriate rescalings of parame¬ 
ters [T5]. A useful rough estimate is that the scattering 
cross section goes as cr = 100a(/ = 100/|(i?) 

A cosmologically interesting level of self-scattering re¬ 
quires cr/rriH ~ l.lb/GeV = 2800 GeV-^ [47] in order 
to address the structure formation problems of cold dark 
matter. This corresponds to a gauge coupling of 

ag = 0.2 /2(i?)(mH/GeV)-3/2 (52) 


The criterion (52) can be compatible with the ioniza¬ 
tion constraint (19) if tor is sufficiently small, 


mil ^ 14 GeV 



(53) 


obtained from eliminating ag from the two relations. 
Very large values of R would be unnatural in our model, 
since it would require a fine-tuned cancellation between 
two contributions to mf = (m^ — yia)^ -I- (j/ 20 ')^, as well 
as a small value of j/ 2 - An accidental cancellation at the 
level of i? = 10 would allow for tor as large as 28 GeV. 


5.2.1. H 2 self-interactions 

In the H 2 phase of the theory, the size of the atom is 
determined by the length scale a that describes the vector 


= 0.19i?2/4(i?2)(^) 


- 3/2 


(54) 


When combined with the constraint (51) on the ioniza¬ 
tion fraction, the result is 


ruR, <6.9GeVi?2/4/3■'^^ 


(55) 


which is similar to the expression found for the H case. 
The primary difference here is that large values of R 2 
can be obtained without fine-tuning of model parame¬ 
ters, allowing for a larger natural range of masses consis¬ 
tent with both the ionization fraction and self-interaction 
constraints. (Notice that —>■ 1 as i ?2 becomes large.) 

Even with a moderate hierarchy R 2 = 10, we can reach 
masses as large as tor^ ~ 70 GeV. 


5.3. Laboratory millicharge searches 

Pair production of is possible in accelerator ex¬ 

periments from the coupling of the photon to the dark 
matter millicharge. The resulting constraints on e are 
quite weak in the mass range relevant for our model, 
mu ^ 1 — 100 GeV, as we show in fig. The existing 
constraints are taken from tables in ref. [IS] for the ASP 
and trident production limits, the E613 beam dump limit 
[49] . ALEPH limits on the Z decay width [50] and a re¬ 
cent GMS search for particles of charge 1/3 or 2/3 [5T] . 
We also show the reach of a new proposed experiment 
for LHC (dashed curve) [52]. These constraints are con¬ 
siderably weaker than that coming from direct detection, 
fig- (lU), which is replotted as the dashed curve on fig. 

On y at low (tor < 4 GeV) or high (tor < 100 GeV) 
masses, outside the sensitivity of direct detection, do they 
become dominant. 

Possibly more signihcant constraints on millicharged 
particles arise from searches for exotic isotopes, bound 
states of normal nuclei with the charged DM constituents. 
Very stringent limits on the concentration of heavy iso¬ 
topes of hydrogen or oxygen from sea water have been 
derived; for example ref. [53j l obtains an upper bound 
of 10“^® for the concentration of anomalously heavy H. 
These experiments assume integer-charged ions, but a 
recent experiment geared toward millicharged particles 
with e > 10“3 set a limit of 10“^^ on the abundance per 
nucleon. Naively such results would seem to rule out al¬ 
most any values of e > 10“® such that the binding energy 
Et = i(ae)^TOp (for anomalous H) exceeds kT at room 
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FIG. 7: Solid curves: existing collider constraints on mil- 
licharge versus mass; dashed curve: expected reach of experi¬ 
ment proposed in ref. [52] . Dotted curves: our direct detection 
limits from fig. depending on choice of Ug = 0.06 (upper 
curve) or Ug = Oion (lower). 

temperature, since we expect some fraction of '0 parti¬ 
cles to remain ionized and thus be able to contaminate 
normal matter. 

However to translate these limits on abundances into 
bounds on e requires many considerations, including the 
expected flux of tp particles, their capture cross section 
on the elements in question, the shielding of the earth 
and the galaxy from charged particles by magnetic fields, 
expulsion of charged particles by supernova winds, the 
process of purification of the samples studied, and the 
question of whether they apply to noninteger charged 
isotopes m- A recent study of these issues was presented 
in ref. [S3]. Here we take the view that there may be 
room for evading the anomalous isotope searches, but 
this question should be revisited if positive evidence for 
millicharges is found. 

6. CONCLUSION 

In this work we have tried to strike a balance between 
simplicity and realism in the construction of an atomic 
dark matter model. Our nonabelian construction is suf¬ 
ficiently rich to explain a unified origin of the massless 
dark photon and charged (under the hidden U(l);i inter¬ 
action) DM constituents Iki as a consequence of symme¬ 
try breaking SU(2)/i —>• U(l)/i by a scalar triplet VEV 
in the dark sector. With the addition of a dark Higgs 
doublet, we have the necessary ingredients to explain the 
asymmetry through leptogenesis, simultaneously with 
the baryon asymmetry. Electric millicharges of '1'^, while 
not a necessary ingredient, can arise naturally through 


heavy states carrying both electric and 11(1)^, charge. 
Higgs portal interactions are also optional, but are al¬ 
lowed by a dimension-4 interaction of 'hi with the dark 
Higgs triplet and its mixing with the SM Higgs. 

The model is mainly testable by direct detection. For 
sufficiently light or heavy constituents, the DM could 
also be discovered in an experiment proposed for LHC to 
probe millicharged particles. It can accommodate strong 
DM self-interactions as suggested by problems of ACDM 
simulations to correctly predict the small-scale structure 
of galaxies, if the dark atoms are not too heavy. Because 
of the requirement ag > Oion, needed to make the ion¬ 
ization fraction in the dark sector sufficiently small, the 
symmetric component of the dark matter is highly sup¬ 
pressed due to annihilations into dark photons, making 
any indirect signals too weak to be detected. 

Our model has a number of features that distinguish 
it from simplified atomic dark matter models. For exam¬ 
ple in the latter, the ratio R of the masses of the atomic 
constituents (which plays an important role) can be arbi¬ 
trarily large, whereas here it is naturally of order 1, and 
requires fine-tuning to be much greater. 

If the new Yukawa coupling yi exceeds the gauge cou¬ 
pling g, the stable dark matter particles can be the lighter 
fermion Iki and the doubly charged (under U(l)/j) vec¬ 
tor boson B , leading to novel three-body bound 

states, where the mass ratio of the constituents mi/niB 
could be large without tuning of parameters (other than 
the usual hierarchy problem of light bosons). The prop¬ 
erties of these unusual atoms for direct detection, as well 
as for DM self-interactions, are qualitatively similar to 
those of the more conventional two-constituent atoms. 
This demonstrates a loophole for strong neutron star con¬ 
straints on asymmetric bosonic dark matter, since the 
dark Coulomb repulsion prevents Bose condensation in 
this model. 

For future work, these models suggest a potential novel 
signal for direct detection, due to the possible simultane¬ 
ous presence of both dark atoms and a subdominant com¬ 
ponent of ionized or symmetric constituents. This would 
allow for the detection of both types of dark matter, typ¬ 
ically having similar but distinct masses and interaction 
cross sections. Our analysis of leptogenesis as a common 
origin of the visible and hidden asymmetries is approx¬ 
imate, and it might also be interesting to undertake a 
more refined treatment for future studies. 
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